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Abstract 



We present a fermion model that is, as we suggest, a natural 2D analogue of the 
Luttinger model. We derive this model as a partial continuum limit of a 2D spin- 
less lattice fermion system with local interactions and away from half filling. In this 
derivation, we use certain approximations that we motivate by physical arguments. 
We also present mathematical results that allow an exact treatment of parts of the de- 
grees of freedom of this model by bosonization, and we propose to treat the remaining 
degrees of freedom by mean field theory. 
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1 Introduction 



Lattice fermion systems in one dimension can be studied successfully by various different 
methods and are by-now well understood. One powerful approach is to perform a particular 
continuum limit leading to a low-energy effective model which can be solved analytically. 
This limit amounts to linearizing the ID tight-binding band relation at the locations of the 
non-interacting Fermi surface (consisting of two points) and then removing the ultra-violet 
(UV) cutoff. In the simplest case of spinless lattice fermions with short-range interactions 
away from half-filling, one thus obtains the Luttinger model [T] which can be solved exactly 
using bosonization [2] ; see [31 HI E] for closely related pioneering work. It is worth stressing 
that exact solubility means a lot in this case: not only the partition function but also all 
Green's functions of the model can be computed by analytical methods; see e.g. [6] and 
references therein. This method can be generalized to ID Hubbard type systems and is 
the basis of a paradigm for ID interacting fermion systems [7]; see e.g. [8] for a textbook 
presentation. We note that bosonization is based on precise mathematical results; see e.g. 
[9] or pig. 

In this paper, we propose a similar approach in two dimensions. Starting from a 2D 
analogue of the spinless lattice fermion system mentioned above, we propose a particular 
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partial continuum limit that makes this system amenable to an analytical, non-perturbative 
treatment in a finite doping regime away from half filling. This limit leads to a model that, 
as we suggest, is a natural 2D analogue of the Luttinger model. Different from ID, only 
parts of the fermion degrees of freedom of this model can be bosonized and thus treated 
exactly. We propose to treat the remaining degrees of freedom using mean field theory. 
We also argue that there is a finite doping regime away from half filling at which these 
remaining degrees of freedom have an energy gap and that in this regime an exactly solvable 
truncation of the model can be used. Our approach is applicable beyond weak coupling. 
Previous work on bosonization in 2D (TTJ [121 E2 HH EEl HHJ EE3 HE] is shortly discussed at 
the end of Section [7J 

Our derivation of this 2D analogue of the Luttinger model relies on approximations 
that we justify by physical arguments; see Section Our proposal that this model gives 
a low energy effective description of 2D lattice fermions is therefore, from a mathematical 
point of view, a conjecture. At any rate, this model is mathematically well-defined and 
can be treated rigorously. To our opinion, it also has a certain mathematical beauty and 
naturalness; see ( )3~T]) - (l3"2l . The details of our approach are quite involved [HI [20j [21] , and 
the aim of this letter is to concisely present the main ideas and results. 

We define our notation and the lattice fermion system that we take as starting point in 
Section [2j Section [3] gives a summary of our results, including a self-contained definition of 
the 2D analogue of the Luttinger model. In Section 0} we discuss physical arguments and 
experiments that motivate and guide our approach. Our derivation of the model and its 
treatment by bosonization are outlined in Sections and El respectively. Section [7J contains 
final remarks. 



2 Two dimensional t-V model 
2.1 Notation 

We consider a square lattice with M 3> 1 sites and denote fermion momenta as k, k'. We 
write k = (ki,k 2 ) = k + e + + A;_e_ with e± = (1, ±1) / \/2, i.e. k± = (ki ± k 2 )/V2, and 
similarly for other 2D vectors p etc. We introduce a lattice constant a > so that — it /a < 
k\,2 < Tc/a, and our large distance cut-off L (system size) is such that k± G (27r/L)(Z + 1/2) 
and L/a e 4 V / 2(N+ 1 /2), i.e. AT = (L/a) 2 . The set of all such k is denoted as BZ (Brillouin 
zone). We use the symbol p for differences of fermion momenta, and A* is the set of all p 
such that p± G (27r/L)Z. Fermion operators ip(k) are defined for k G BZ and normalized 
such that {?/>(k), ?/^(k')} = 5k,k'[-^/(2vr)] 2 , and we use the same normalization for other 
fermion operators ?/v iS (k) introduced in Section [3) We write ^^(k) short for ^(^^(k) 
etc. We use the abbreviations a = 2\/2a. The symbol 5 a denotes the lattice periodic 
Kronecker delta, i.e. S a (p) = 1 if p G (27r/a)Z 2 and otherwise. Finally, [k + k'] is the sum 
of two momenta in BZ modulo (27r/a)Z 2 , i.e. it is equal to k + k' + (2n/a)n G BZ with 
appropriate n G Z 2 . 
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2.2 Definitions 



The lattice fermion model we consider is the so-called 2D t-V model. It describes spinless 
fermions on a 2D square lattice with Af sites and hopping and repulsive density- density 
interactions between nearest neighbor sites. It is defined by the Hamiltonian 

H tv = H - f iN + H ink (1) 

with the free part 



#o=(^) 2 $>(k)^ (k) (2) 

keBZ 

and e(k) = — 2t[cos(aA; 1 ) + cos(a/c 2 )] the tight binding band relation, 

^=(¥) 2 E^( k ) ( 3 ) 

keBZ 

the fermion number operator, and \x G K the chemical potential. The interaction is 

#int = (tT E ^ t (ki)^(k 2 )^t(k 3 )v;(k 4 )M(k 1 - k,)^ - k 2 + k 3 - k 4 ) (4) 

kj GBZ 

with -u(p) = a 2 V[cos(api) + cos(ap 2 )]/(87r 2 ) the Fourier transform of a nearest neighbor 
interaction. The model parameters t (hopping constant) and V (coupling strength) both 
are positive. The filling parameter is defined as 

where (•) denotes the ground state expectation value. It is in the range < v < 1 with half 
filling corresponding to v — 1/2, and v — 1/2 is referred to as doping. 

The Hamiltonian H t v is invariant under the particle-hole transformation 

V>(k) ->■ -0 f (k + Q/a), z/ — >• 1 — i/ (6) 

for Q = (7r, 7r), up to an irrelevant additive constant. We therefore can restrict ourselves to 
v > 1/2. 



3 Summary of results 

The 2D analogue of the Luttinger model describes six flavors of fermions ^v jS (k) with flavor 
indices r — ± and s — 0, ± and momenta k in different Fourier spaces A* as follows, 

K = {k\k ± e 2 f(z + ±), -Z<k ± <l\ 

l L 1 a a) 

A* =± = |k | k± e ^(Z + |), -oo < k s < oo, 



71 

■- < k- s < 
a 



-}■ 

a J 



(7) 
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We refer to the fermions with s = and s = ± as nodal and antinodal, respectively. The 
Hamiltonian defining this model is 

H = H n + H a (8) 

with the nodal part 

H n =E((x) 2 E E™^ ■i>tjr^)- + 



s=± keA* r=± 

peA* r =± r,r'=± 

and the antinodal part (including nodal-antinodal interactions) 

H a = (x) 2 E J2(~ rCFk + k - ~ ^ : $Ao(k): + 

keA*, r=± 

(l) 2 E ^ 3 E ^,o^-r,o(p) + 9a E X(P)4,0^',,(P 
peA* r =± r,r',s=± 



(9) 



(10) 



(the model parameters t^, cp, g>j, and fi a are specified further below). The J r>s are Fourier 
transformed and normal ordered fermion densities defined as follows, 



4o(p) = (f) 2 E :$,o(k-p)i, Q (k): 
keAS 

4,=±(p) = (f ) 2 E E ^J,.(k-p)^(k + 27me_./S): 



, 2 

keA* nG2 



'11} 



and J} s J r , 8 /(p) is short for J rs (— p)J r / jS /(p). We also use the cutoff functions 

x(p) = 1 if — vr/a < p± < 7r/a and otherwise. (12) 

The colons indicate normal ordering with respect to a non-interacting groundstate Q defined 
by the conditions 

A. j± (k)fi = $ )± (-k)fi = if rfc± > (13) 

and 4>l (k)Q = or ^ r0 (k)fi = if rk + k^ > or < 0, respectively. In the theorem below 
we set J r ,s=±(p) = for \p- s \ > it /a. 

We first outline how the model above can be derived from the 2D t-V model using 
certain approximations (Section [3} . This allows us to determine the model parameters as 
follows, 

c F = 2ta 2 , v F = 2V2tasm(Q), Q = tiv (14) 
0i = 2g 2 = 2V sin 2 (Q)a 2 , g 3 = g 4 = 2Va 2 (15) 
V\ ^ , TA 2( ^, ( ^ 2Q 



fi a = - (4t + -j cos(Q) + V cos 2 (Q) ^1 - ^ (16) 
with < \v — 1/2 1 < 1/4. The key result to bosonize the model above is the following. 



Theorem: (a) The nodal density operators obey the relations 

[J r ,s{p), J r ',s'(p')] = 8r,r'Ss,s'—r^Sp,-p> (^) Vr, S = ± (17) 

and Jr.sipV = Jr,s{—p)- Moreover, J r ^ s (p)Q = Vp such that rp s > 0. 
(b) The nodal Hamiltonian in is identical with 

H n = (l) 2 ^2 (nWF : Jr,sJrAP) '■ +X(p) [fl , l«A-,«^-r,a(p) + #2 ^ J r , s J r ' _ s (p)] J (18) 

pe A* r,s=± r'=± 

and well-defined provided that V < 47rt/ sm(Q). 
(Proof outlined in Section |6j) 

As discussed in Section El his theorem provides the means to exactly map H n to a 
bosons Hamiltonian. Thus, the model above is equivalent to a model of non-interacting 
bosons coupled linearly to the antinodal fermions. It is therefore possible to integrate out 
the nodal fermions exactly and thus obtain an effective antinodal model. We also discuss 
how these results can be used to obtain physical information about the 2D t-V model. 



4 Physics motivation 

Our aim is to rewrite and modify the 2D t-V model such that, (i) the resulting model can 
be treated by bosonization, (ii) the low energy physics is changed as little as possible. For 
that we adopt the following hypothesis that has been successful in ID: 

HI: There exists some underlying Fermi surface dominating the low energy physics, and 
we can modify, ignore or add degrees of freedom far away from this surface (in the latter 
two cases we need to correct the definition of doping, of course). 

We note that, at half filling and without interactions, the relevant Fermi surface is the 
square \ki ± k%\ = ir/a (the large diamond in Figure 1). We assume that the underlying 
Fermi surface remains close to this square even in the presence of interactions and away 
from half filling. We observe that the band relation has a qualitatively different behavior 
in different regions of the Brillouin zone close to such a surface. To make this explicit, we 
select six representative points Q r< s/ a labeled by two indices r = ± and s = 0, ± as follows, 

Q+,o = (7r,0), Q_, = (0,7t) 
Qr s = (rQ, rsQ) for r, s = ± 

for some Q ~ 7r/2 (the dots in Figure 1). We then approximate the band relation close to 
these points by truncating the Taylor series at the lowest non-trivial order, e(Q ryS /a + k) w 
e r>s (k) with 

e r> o(k) = —rcpk + k_, e n ±(k) = — 4tcos(Q) + rvpk± (20) 

and the constants Cp and vp in (|14|) . At half-filling we expect Q = 7r/2, and, for symmetry 
reasons, the underlying Fermi surface should always contain four such points Q n ±, r = ±, 
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for some value of Q. Thus, the band relation is hyperbolic close to the vertices Q r ,o/a of 
the square surface, and it is linear with a constant Fermi velocity Vp close to the midpoints 
of the sides Q r ,±/a. We use the terminology of experimental physicists studying cuprate 
superconductors with ARPES [23] and refer to the regions in Fourier space close to Q r fi/a 
and Q r ,±/a as antinodal and nodal, respectively. 

As will become clear in the next section, Hypothesis HI is not enough to justify the 
approximations we make. For that the following somewhat stronger hypothesis is needed. 

HI': The low energy properties of the model are not (much) changed if we modify the band 
structure and the interactions for fermion degrees of freedom far away from the six Fermi 
surface points Q rjS /a in / fiPj) . 

There are physical arguments suggesting that, at half filling and sufficiently large cou- 
pling, the 2D t-V model has a charge density wave (CDW) groundstate which is insulating; 
see e.g. Section IV. A in Ref. |22j. In mean fields theory this state is characterized by a 
CDW gap A > changing the band relations to ±A/e(k) 2 + A 2 , and this corresponds to 
a fully gapped underlying Fermi surface. Moreover, ARPES results on cuprate supercon- 
ductors show that there exists an interesting doping regime away from half filling where, in 
these materials, the antinodal fermions are gapped while the nodal fermions have no gap 
but Fermi surface arcs [23]. This suggests that the nodal- and antinodal fermion degrees 
of freedom can have different physical behavior, and it motivates us to rewrite the 2D t-V 
model so that nodal- and antinodal fermion degrees of freedom can be easier treated by 
different computation methods. 

The model thus obtained can be further simplified by the following hypothesis motivated 
by the ARPES results mentioned above. 

H2: There exists a finite doping region away from half filling where the nodal points on the 
underlying Fermi surface move to Q r ,±/ a with Q > tt/2, and in this regime the CDW gap 
is absent in the nodal regions while it is still present in the antinodal regions. 

In this regime, the parameter Q is determined by doping. This hypothesis is also 
suggested to us by results from renormalization group studies [2H [25] and mean field theory 
[261 [27]. 

5 Derivation of the model 

We now outline a derivation of the model in fl51)- f[T6l) from the 2D t-V model, emphasizing 
the approximations that are made. 

We divide the Brillouin zone in six regions as shown in Figure 1 (the areas enclosed by 
boldface lines), and we label the fermion degrees of freedom in these regions as follows, 

4,s(k) =^([Q r , s /a + k]) (21) 

where the momenta k are now restricted to the small regions A* s containing only momenta 
satisfying — n/a < k± < n/d for s = 0, and — tt/cl < k_ s < n/d, —n/(2a) < ki + r{Q — 



6 




Figure 1: Division of the Brillouin zone in six regions A* s , r = ±, s = 0, ±. The six dots 
mark the points Q r , a / a in ffT9l for filling v = 0.55. The regions A* s are equal to the six-sided 
polygons (for s — ±) and small diamonds (for s = 0) shifted such that the corresponding 
points Qr^/a coincide with the origin of the coordinate system. Shown are also the half- 
filled square surface (large diamond) and the assumed interacting Fermi surface for the 
nodal fermions (dashed lines). 

n/2)/a < n/(2a), -7r/(2a) < k 2 + rs(Q - n/2)/a < 7r/(2a) for s = ±; see Figure 1. We also 
assume 7r/4 < Q < 37r/4 and Q G \/27r(a/L)M (the latter restriction becomes irrelevant in 
the limit L/a — > oo). We thus trade different regions in the Brillouin zone for six different 
fermions flavors r, s. Note that A?j_ = Aq in ([7]), and the regions are defined such that 

E/( k ) = E E /([Qr,./a + k]) (22) 

keBZ r,s kGA*. >s 

for any function / on BZ where, here and in the following, sums over r and s are over r = ± 
and s = 0, ± unless stated otherwise. In particular, = J2rs-^r, s with 

N r,s= (t) 2 E (23) 

keA* s 

and ifo = X] r s(2 7r /-^) 2 SkeA* e (Qr,s/^ + k)V ; J s '?/ ; rs (k). We now make our first approxima- 
tion: 
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Al: Taylor expand the exact band relations in the vicinity of the Fermi surface points 
Q, r ,s/ a °> n d keep only the leading non-trivial terms, i.e. replace e(Q r>s /a + k) in H above 
by'e r>s (k) in (H2|). 

We thus obtain 

H o = J2(t) 2 E (24) 

r,s keA* 

In this approximation, we neglect terms of higher order in ak. We expect that the low 
energy properties of the model are dominated by states close to the Fermi surface points 
for which these higher order terms are negligible. 

To modify the interaction part of the Hamiltonian in a way suitable for our purposes 
we use ( 122]) to rewrite 

#int = {tY E E i t 1 , Sl (ki)^ 2 , S2 (k 2 )^ 3jS3 (k 3 )i 4iS4 (k 4 )«(Q riiSl - Q r2iS2 + kx - k 2 ) 

r,s kj-eA* >s 

x<5 a (Q riiSl Q r2jS2 + Qr3,s3 Qr4,s4 + ki k 2 + k 3 k 4 ) 

(25) 

and make our second approximation: 

A2: Replace the interaction Hamiltonian in / T^3)) by 

H **= (x) 6 E E ^ 1 ,, 1 (ki)^ 2 , S2 (k 2 )^3, S3 (k3)^ 4iS4 (k 4 )«(Q riiSl -Q r2 , S2 ) 

(26) 

^^a(Qri,si Qr2,S2 4" Q»"3,S3 Q^^^ki— k 2 +k3— k4,0- 



Note that Q riiSl Q r2iS2 -|- Q r . 3jS3 Qr4 

iS4 is either or a "large" vector whose length is 
of order n/a. We expect that the most important processes for the low energy properties 
of the model are those where the kf are so small that the Kronecker delta in (1251) can be 
non-zero only if k x — k 2 + k 3 — k 4 = 0, and these processes are not much affected by A2. 

To simplify (1261) we have to find all solutions of <5 a (Q riiSl — Q r2iS2 + Qr 3 ,s 3 — Qr 4 ,s 4 ) = 1- 
We assume Q ^ it/2 and obtain after straightforward computations [19] 



H ^ = (l) 2 E [93^2 PlfiP-rfiiP) + 9l E Pl,sP-r,siv) 



2 2^l9s 



+92 Yl Pt,sPr',-s(P)+94 E Pt,0Pr>, S (P 
r,r',s=± r,r',s=± 



(27) 



with 



(P) = (t) 2 E ^(ki)VUk 2 )<Wk 2 (2f 



k,-eAj 



and the coupling parameters in (fl5l) . 



The operators in have the natural physical interpretation as Fourier transformed 
density operators. It is remarkable that we only get interaction terms of density-density 



8 



form. Our assumption Q ^ 7r/2 is important since otherwise one obtains additional terms 
that cannot be treated in a simple manner by bosonization, similarly as in ID [8]. Such 
terms might cause a nodal fermion gap. This is consistent with our Hypothesis H2 in 
Section H] since Q = tt/2 corresponds to half filling, as shown below. 

To obtain a model that can be bosonized we want to remove the UV cutoff for the 
nodal fermions orthogonal to the Fermi surface arcs. Before we can do that we have to 
normal order all fermion operators with respect to some suitable reference state Q ( "Dirac 
sea") with fixed fermion density v ^ 1/2. We choose this state such that the nodal fermion 
states are filled up to the Fermi surface through the points Q r ,±/ a parallel to the half- 
filled one (the four dashed lines in Figure 1). For simplicity we also assume that the 
antinodal fermions are half-filled (one can show that this follows from our Hypothesis H2), 
but this is not essential (see [20] for a more general treatment). This reference state can be 
defined by the conditions in ([T31 ff. By simple geometric considerations we find that the 
contributions of the different regions r, s to the total filling v of this state are v r $ = 1/16 
and v r ,s=± = (Q/n — 1/8) /4, respectively, and thus v = J2 rs u r,s — Q/tt. 

Normal ordering of the Hamiltonian H — fj,N + H int amounts to the following: rewrite 
H mt in terms of the normal ordered fermion densities 

J r ,s(p) = ■■ Pr,s(P) '■ ( 29 ) 

and normal order H and N, dropping an irrelevant additive constant. An important 
consequence of this and Approximation Al is a change of the chemical potential term — /j,N 
in the Hamiltonian to — J2rs^r,s '■ N r>s :, and we find by straightforward computations 
the following renormalized chemical potentials: /i rj0 = // — (2Q/tt)V = fi a and /i r)S= ± = 
H + At cos(Q) — [2Q sin 2 (Q)/ it + cos 2 (Q) +cos(Q)/4]V = // n . To have an underlying Fermi 
surface as assumed in H2 we set /i„, = 0. This fixes \i and implies (TTBI) . Note that \i a changes 
sign under the transformation Q — > n — Q, as required by invariance under particle-hole 
transformations. 

By assumption the groundstate expectation values of : N r>8 : all are zero in the parameter 
regime of interest to us, and thus the filling factor of the groundstate is identical with 
v — Q/tt. This fixes Q as in (|14p . 

We now can partially remove the UV cutoff for the nodal fermions. We do this by a 
sequence of three approximations: 

A3: Replace in the interaction part of the Hamiltonian the normal ordered nodal densities 
Jr,s=±(p) by x(p) Jr,s=±(p) with the cutoff functions x i n (HIP- 

A4: Replace the nodal Fourier space regions A* s=± by the sets A* =± in |?|), i.e., drop the 
restriction on k s . 

A5: Replace in the interaction part of the Hamiltonian the normal ordered nodal fermion 
densities 

l,s=±(p) = (f ) 2 £ :^(k-p)i, s (k): (30) 
keA* 

by the ones in / TO]) , i.e. add the umklapp terms n^O parallel to the Fermi surface arc. 
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We thus obtain the model in (!H|)- (|T6|) . 



Note that in approximations A4 and A5 we add terms to the nodal fermion densities, 
and we therefore increase the number of interaction terms in the Hamiltonian. Approx- 
imation A2 partly compensates for this. To be more specific: before A3, the number of 
terms included in the nodal fermion density J r ,s=±{p) decreases linearly with \p + \ and |p_| 
until \p±\ = 2ir/a when it becomes zero. However, after A4 and A5 this number of terms 
becomes infinite and independent of p. It therefore is natural to restrict the interaction 
terms involving nodal fermions densities by imposing a cutoff on p. We do not know a 
physical argument that fixes this cutoff in a unique manner. However, we believe that the 
properties of the model are not sensitive to changes of this cutoff. We therefore choose here 
a cutoff that is particularly simple; see [19] for a more general treatment. 

Our arguments above can be generalized to lattice models with more general band re- 
lations, and our choice for the regions A* can be generalized by introducing an additional 
parameter allowing to change the widths of the nodal regions A^_ [19]. Moreover, it is not 
necessary to use the simplified band relations e r> o and interaction vertices for the antin- 
odal fermions, and we only do this here for aesthetic reasons. It also is possible to relax 
Hypothesis H2 [20J, as already mentioned. 

It would be interesting to check in more detail if the approximations above indeed do 
not change the low energy properties of the model. This can be done, in principle, using our 
results: one can bosonize also the terms that we added or dropped in our approximations, 
and one should be able investigate the effect of these terms on the groundstate using 
renormalization group methods. This is an interesting project for the future. 



6 Bosonization and partial exact solution 

The arguments in this section can be made mathematically precise |19[ l2~Tj . 

We first outline a proof of the theorem in Section [3j For that it is useful to formally 
write the model in (|9"|)- (|TT|) in position space as follows (we suppress arguments x below), 

H n = / d 2 X ^ '■ rV F 1pt,a {-ida)i>r,8 +^1^ Jr, S J-r,s + 92^ J r, S Jr'-sj ■ (31) 

s=± r r r,r' 

for the nodal part (d± means d/dx±), and 

H a = d 2 X : 4>lfi( rC F d + d - ~ ^a)A,0 + 93^2 J r,Q J -r,0 + 94, J rfl J r',s) : ( 32 ) 

r r s=±,r,r' 

for the antinodal including the mixed parts, with J rjS (x) =: ■0j a (x)^' ra (x) : the normal 
ordered fermion densities. The proper interpretation of these formulas follows from the 
precise definition of the model in Fourier space: different fermion flavors ip TyS come with 
different Fourier space regions A*, and thus the spatial variable x in 

W, s (x) = ( 2 f f #rM ex P( ik ' x ) (33) 

keAj 
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lives on different spaces: for the fermions ?/v,+ (x) only the variable x + is continuous while 
x_ lives on a ID lattice with lattice constant a and length L, and similarly for ?/v _(x) with 
x + and X- interchanged (note that the umklapp terms with n 7^ in (TIT]) are needed for 
this interpretation to be consistent with Fourier transformation). Thus the integrals here 
should be (partly) interpreted as Riemann sums [19]. In particle physics parlance, 1/5 is a 
UV cutoff needed to give a precise mathematical meaning to the model. 

Thus ?/v j+ (x) in ( I3T1) can be treated as a collection of ID Dirac fermions where x + plays 
the role of the ID space coordinate and x_ of a flavor index. We therefore can use the 
standard mathematical results of ID bosonization; see e.g. [9] or [10J. These imply that the 
densities J rj+ obey the following commutator relations, 

[J r ,+(x), Jr',+(y)] = S r yr— -d+S 2 (x - y) (34) 
where 5 2 (x — y) here means 5(x + — y+)5 x _ ty _/a, and 

r J d 2 x :V4, + (-id+)VV,+ : = a J d 2 x :J 2 + : (35) 

where f d 2 x stands for f dx + ^2 X _ ®- These and similar formulas for ip r _ imply (TTT|) and 
( TT8"j) . The other relations stated in part (a) of the theorem in Section [3] follow from well- 
known corresponding results in ID. 

We now describe how to exactly diagonalize the bosonized nodal Hamiltonian H n . For 
simplicity we continue our discussion in position space and ignore UV cutoffs that do not 
affect the results we mention (a more detailed solution, including all cutoffs explicitly, can 
be found in [19J). The commutator relations in (|34"|) imply that the operators 

9±$± = Vrii(J+,± + J_,±) 
U± = y/nd,(-J + ^ ± + J_ ±) 

obey the commutator relations of standard 2D boson fields: [$±(x), n±(y)] = i5 2 (x — y) 
etc., and by straightforward computations one finds that the nodal Hamiltonian can be 
expressed in terms of these bosons as follows, 



H n = ^jd 2 x: (g[(l - 7 )n 2 + (1 + 7 )( W 2 ] + 2 7 (d + $ + )(«9_$_; 



(37) 



where 7 = V sm(Q) / (Ant) . It is important to note that this Hamiltonian is positive definite 
only if 7 < 1, and this implies an upper bound on the interaction strength mentioned in the 
theorem in Section [31 similarly as in ID. Fortunately, this bound is satisfied for interesting 
parameter values (V/t from up to 12 or so). 

We thus obtain an exact representation of the 2D analogue of the Luttinger model as 
a system of non-interacting bosons coupled linearly to the interacting antinodal fermions. 
The boson Hamiltonian in ( 1371) can be diagonalized by standard methods, and we find 

H n = £ n + ( 2 f) 2 E w.(p)6l(p)6.(p) (38) 

p s=± 
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with some computable constant £ n , the dispersion relations 

w±(p) = |P| 2 ± \/|p| 4 - (1 - [7/(1 + l)] 2 )(2 P+ p!y, (39) 

and standard boson operators b s (p), i.e [6 ± (p), 6±(p')] — <5 2 (p — p') etc. Note that the 
anisotropy of the non- interacting systems is reduced by 7 7^ 0. Using that it is straightfor- 
ward to compute the exact energy eigenstates, energy eigenvalues, and partition function 
for the model defined by the Hamiltonian in (Q. 

One can also exactly integrate out the bosons and thus derive an effective action for the 
antinodal fermions. We found that, in a good approximation, this action can be represented 
by an effective Hamiltonian that has the same form as H a in (flOl) but with the coupling 
constants changed to [T5J 

g 3 = 2Va 2 (l- . V , V s 4 = 0. (40) 

\ 2 sm(Q)|27rt + V sm(Q)\ J 

It is possible to compute the CDW gap A oc (ijr + ip_ ) for this latter effective antinodal 
model using mean field theory [20]. We found that (:N a :) = (:N n :) =0 and A > 
are indeed independent of the parameter /j, a for intermediate values of V/t and in some 
Q-interval around Q = 7r/2 [20], in agreement with Hypothesis H2 in Section H]. 

We believe that, after bosonization, it is possible to make mathematical sense of the 
continuum limit a — > of the nodal Hamiltonian, but we expect that this requires non- 
trivial additive and multiplicative renormalizations, similarly as in ID [28]. This would give 
further support to our conjecture that the approximations in Section \5\ do not affect the 
low energy properties of the model. 



7 Final remarks 

We argued that there exists a parameter regime away from half filling where the antinodal 
fermions Vv,o are gapped. We propose that, in this regime, the nodal Hamiltonian in (Q 
alone accounts for the low energy physics of the 2D t-V model. This nodal model is exactly 
solvable by bosonization and, in particular, it is possible to compute all its Green's function 
by analytical methods. It would be interesting to do this and thus confront this model with 
experimental results on 2D correlated fermion systems. We hope to come back to this in 
the near future [21]. 

While the 2D analogue of the Luttinger model becomes particular simple in the regime 
where the antinodal fermions are gapped, it can be used also for other parameter values 
[20J. We believe that a mean field treatment of the antinodal fermions is appropriate to 
compute phase diagrams. However, in parameter regions where mean field theory predicts 
that the antinodal fermions are gapless, a treatment of the antinodal Hamiltonian in (fTUj) 
beyond mean field theory would be desirable. 

Our approach can be straightforwardly generalized to 2D lattice fermion models with 
more complicated band relations [T9~| |2"U] and spin [2T], similarly as in ID [8]. 
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It is known that the truncated nodal model in (jUj) with g 2 = has "Luttinger liquid" 
behavior after [29J but not before Approximation A5 [30]. However, we believe that our 
model with g 2 > is less sensitive to Approximation A5 than the truncated model with 
g 2 = 0. The reason is that the boson propagator has better decaying properties in Fourier 
space for g 2 > than for g 2 = (since the boson dispersion relations in (I3"9~j) are quite 
isotropic, whereas for g 2 = they are cu±(p) oc \p±\ independent of p T ). 

Although we suggest that the model in ()8l)-(ITTT) is a 2D analogue of the Luttinger model, 
we emphasize that we do not state that this model has "Luttinger-liquid" behavior [15] : if 
or not this is the case is a delicate question and remains to be seen. We plan to come back 
to this in the near future |21j . 

We finally discuss previous work related to ours. The pioneer in higher dimensional 
bosonization is Luther [11] . To our knowledge it was Mattis [12] who first made the im- 
portant observation that a 2D model similar to our nodal fermion model (i.e. no antinodal 
fermions) is exactly solvable by bosonization; see also Ref . [T3| H3] for previous work on Mat- 
tis' model. Our work can be regarded as a particular implementation of the idea promoted 
by Anderson [T5] that 2D interacting fermions systems can be understood by bosonization 
in each direction of some Fermi surface. An earlier implementation of this idea and starting 
with a square Fermi surface is by Luther [16], but different from him we bosonize chains 
in position space (rather than Fourier space) and thus can avoid certain approximations in 
the treatment of the interactions. Various other implementations of the idea to bosonize 
2D fermion systems appeared in the literature before but seem to differ in detail from ours; 
see e.g. [T7J [T5] and references therein. We finally mention other work that was impor- 
tant as inspiration for us, namely renormalization group studies of the 2D Hubbard model 
|24[ [25] and work by Schulz [26] emphasizing the significance of the antinodal fermions for 
2D lattice fermion systems. 
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